First-principles Simulations of the stretching and final breaking of Al nanowires: 
Mechanical properties and electrical conductance 
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The evolution of the structure and conductance of an Al nanowire subject to a tensile stress 
has been studied by first-principles total-energy simulations. Our calculations show the correlation 
between discontinuous changes in the force (associated to changes in the bonding structure of the 
nanowire) and abrupt modifications of the conductance as the nanowire develops a thinner neck, in 
agreement with the experiments. We reproduce the characteristic increase of the conductance in the 
last plateau, reaching a value close to the conductance quantum Go = 2e 2 /h before the breaking of 
the nanowire. A dimer defines the contact geometry at these last stages, with three channels (one 
dominant) contributing to the conductance. 

PACS numbers: 73.63.b, 62.25.+g, 73.63.Rt, 68.65.-k 



The electrical and mechanical properties of metallic 
nanowires have received a lot of attention Q. Although 
point contacts have been studied for many years, only 
recently the gentle control of the distance between two 
metals using an STM-AFM § or a mechanically con- 
trolled breaking junction (MCB J) || has allowed the ex- 
perimental characterization of atomic contacts and the 
observation of quantum effects in both the conductance 
and the forces In a pioneering work, Scheer et al.Q 
have shown, analyzing the superconducting properties of 
an atomic contact, how the transport properties of the 
system just before the breaking point depend on a few 
channels that they related to the atomic orbitals of the 
contact. 

The formation of necks and atomic contacts in 
stretched metallic nanowires has been analyzed theoreti- 
cally using different approaches. The nanowire deforma- 
tion can be studied by molecular dynamics simulations 
using either classical || [| or effective-medium theory 
potentials [M ||. First principles calculations based on 
Density Functional theory (DFT) §, |o| provide a more 
accurate description of the mechanical properties and the 
electronic structure needed for the calculation of the con- 
ductance, but the large computational demand restricts 
most of the applications to the analysis of model geome- 
tries for the contact (e.g. monoatomic chains). Up to our 
knowledge, the most complete calculation of the deforma- 
tion of a metallic nanowire has been presented by Naka- 
mura et al Q, who analyzed, using DFT calculations, a 
Na nanowire with 39 atoms. In this simulation, the wire 
is elongated in steps of 0.2 or 0.4 A until it reaches the 
breaking point. The conductance is determined using the 
Landauer-Buttiker formula, where the transmission ma- 
trix is calculated from the self-consistent electrostatic po- 
tential using scattering techniques Q. This calculation 
showed how the nanowire deformation was accompanied 



by a rearrangement of the atomic configuration, that in- 
troduces also jumps in the forces and the conductance of 
the system. 

Addressing this complex problem with a fully- 
converged first-principles description would be still too 
computationally demanding. One possibility is to stick 
to accurate plane- wave (PW) DFT methods, carefully re- 
laxing the conditions for convergence jl3| (basis set cut- 
off, k-point sampling, etc). A second alternative is to 
resort to local orbital MD-DFT methods, specially those 
devised with the aim of computational efficiency, that 
allow first-principles studies of much more complex sys- 
tems. The formulation in terms of local orbitals has an 
added value, as the transport properties can be easily cal- 
culated from the resulting (tight- binding or LCAO) elec- 
tronic hamiltonian using non-equilibrium Green's func- 
tion techniques JlJ]. Thus, efficient local-orbital MD- 
DFT methods are probably the best available tools for 
a first-principles analysis of complex nanowires. In this 
work, we have studied the evolution of the structure and 
the conductance of an Al nanowire as a function of the 
stretching by means of a fast local-orbital minimal ba- 
sis MD-DFT technique (Fireball96)Jl|. This technique 
offers a very favorable accuracy/efficiency balance if the 
atomic-like basis set is chosen carefully [|l6|. This ap- 
proach is complemented with PW-DFT calculations (us- 
ing CASTEP ) performed at the critical points of the 
nanowire deformation, as discussed below. Closely re- 
lated approaches jl9|, |2(], |l|, |2^], based on the com- 
bination of ab-initio calculations on local-orbital basis 
and Keldysh-Green function methods, have been recently 
applied to characterize some ideal geometries (mainly, 
atomic chains) for Al and Au nanowires. 

We have studied the system shown in figure [l]A, hav- 
ing 48 atoms -12 of them in the wire-, embedded in 
an Al(lll) surface having a (3x3) periodicity. We also 
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FIG. 1: Ball-and-Stick model of the structure of the Al nanowire for different steps of the stretching process (see fig. g). The 
atoms involved in the important bonding rearrangements related to discontinuous changes in total energy, force and conductance 
are labelled 1 to 5. 



impose periodic conditions in the direction perpendicu- 
lar to the surface, joining artificially the last two layers 
of the system. The stretching of the system is simulated 
increasing the distance between the two limiting layers 
by steps of 0.1 A. After each step, the system is allowed 
to relax toward its configuration of minimum energy; in 
this relaxation, only the atoms located in the last two 
layers remain fixed, meaning that 30 atoms are allowed 
to relax. We have checked the validity of the calculation 
performed with Fireball96, repeating similar calculations 
with CASTEP at points of the deformation where rear- 
rangement of the atoms appeared (see the discussion be- 
low) . Figure |l| shows snapshots of the nanowire geometry 
along the stretching path; different profiles correspond 
to the points labelled A to F in fig. |], where the total 
energy per atom is shown as function of the stretching 
displacement. This figure || displays the points where 
the nanowire suffers an important rearrangement in its 
geometry: these jumps occur between points B and C, 
and between points D and E. Fig [j] shows the initial 
geometry (A), the geometries before and after the first 
(B,C) and second (D,E) jumps, and the structure close 
to the final breaking point (F). It is important to no- 
tice that the main rearrangement of atoms occurs in the 
layer formed by atoms 1,2 and 3. In particular, along 
the path AB, the atom 3 starts to move out of the initial 
1-2-3 layer; at the first jump, this atom takes an indepen- 
dent position linking a lower layer and the two atoms 1-2. 
Along the path CD, the geometry of the point C does not 
change too much, the atom 2 moving slowly toward atom 
3; eventually, at the second jump, atom 2 forms a dimcr 
with atom 3, and the nanowire reaches a geometry very 
similar to the one formed at the breaking point. In par- 
ticular, along the EF path the system evolves basically 
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FIG. 2: Total energy of the Al nanowire as a function of the 
stretching displacement, Ad. 

changing mostly the distance between atoms 2 and 3. 

We have checked these calculations, obtained using 
a Fireball code, by recalculating geometries around the 
jumps (BC and DE) using CASTEP: these calculations 
confirm the validity of the Fireball results; we have 
found, however, some minor differences. For instance, 
in CASTEP we find the system to present the first jump 
for Ac? = 0.9 A, namely, a step before the jump found 
in Fireball; the system evolves, however, across the first 
jump to the same geometry calculated in Fireball. 

Figure [?] shows the force along the nanowire stretching. 
It is very satisfactory that the first jump seems to appear 
for a maximum in the force; along the path CD the force 
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FIG. 3: Normal force along the stretching path. The inset 
shows the force between n.n. atoms in Al crystal defined 
as the derivative of the energy per atom with respect to the 
n.n. distance during hydrostatic strain, divided by 6 bonds 
assumed in the f.c.c. structure. 



FIG. 4: Total differential conductance (in units of the conduc- 
tance quantum) and channel contribution along the stretching 
path. 



is practically constant and along the final deformation, 
path EF, the force increases except for some small fluctu- 
ations, until reaching a maximum at the point where the 
wire breaks. It is interesting to compare the forces shown 
in H with the one calculated at the break point for an Al- 
crystal: this is defined by the maximum derivative of the 
energy per atom with respect to the nearest-neighbor dis- 
tance for an Al-crystal, divided by the number of bonds 
each atom is assumed to form (6 bonds per atom in a 
f.c.c. structure). Our calculations yield 0.7 nN for the 
maximum force, and this quantity should be compared 
with the forces of fig |^ that are typically larger than that 
value: at the nanowire break point the force is around 1.2 
nN, 1.7 times the bulk value. 

We have calculated the electrical conductance of the 
nanowire using a Kcldysh-Grccn function approach based 
on the first-principles tight-binding Hamiltonian ob- 
tained from the Fireball code, at each point of the defor- 
mation path. In this formalism, we rewrite this hamil- 
tonian describing the system as Hi + H 2 + T12, where 
the total system is splitted into the two parts, 1 and 2, 
T12 defining the coupling between both. Typically, we 
use the thinnest part of the nanowire to define the inter- 
face between these two subsystems. Then, the differential 
conductance E3 is given by: 
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[l-T 12 g r 22 (E F )T 21 g r n (E F )Y 
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G - dV 

Aire 



(1) 



— ImTr[T l2 p 22 (E F )D r 22 (E F )T 12 p n (E F )D^(E F )] 

where pn and p 22 are the density matrices associated 
with sides 1 and 2, respectively, while: 



and is given by a similar eqn.; and g 22 are the 
retarded Green function of the decoupled sides, 2 and 1. 

Using eqn. [j], we calculate the differential conductance 
along the stretching deformation (fig. ^From this fig- 
ure we see how this conductance changes dramatically at 
the points B-C and D-E: in both cases the conductance 
jump coincides with the force jump (fig ^), and is deter- 
mined by the atomic rearrangement the system suffers at 
those points. The conductance behavior shows two other 
important features: (i) the nanowire breaks at a point 
where the conductance is very close to the quantum unit, 

G w 0.95 (ii) second, before the breaking point, 

the conductance takes values a little smaller than the 
last one, with a negative slope. These two features are in 
good agreement with the experimental evidence |l], |2j| . 
Morever, at the breaking point the system shows a ge- 
ometry slightly different to the one suggested in 0]. In 
particular, we find that an Al-dimer is responsible for 
the conductance properties the nanowire shows before it 
breaks. 

Finally, we have also addressed the question of how 
many different channels contribute to the conductance. 
This has also been analyzed using eqn [l]; this eqn, can 
be rewritten, using the cyclic property of the trace in the 
form: 



(3) 



where i — p\\D\^T\ 2 p 22 p3. By diagonalizing this trans- 
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fer matrix, we find the channels contributing to G. 

Figure |] also shows how different channels contribute 
to the total differential conductance. Near the break 
point, we find that only three channels yield appreciable 
contributions to the total current, although there appears 
a predominant mode, very much in agreement with the 
results published by Scheer et al |Q. Along the trajec- 
tory CD, before the last jump, the system also presents 
three predominant channels, but, in this case, all the 
three yield important contributions to the current; no- 
tice that in this case, the total conductance is around 2 
(^f-), with a channel contributing practically (%-), and 

two other channels, each one contributing if^f-)- In the 
initial path, AB in fig 1, we find six channels contribut- 
ing with a total conductance of around 3(^-). From 
this analysis we conclude that the channels are not either 
open or closed; their particular contribution depends very 
much on the geometry of the nanowire contact. We can 
only say that channels tend to close along the nanowire 
deformation, and that once a channel is closed by the 
stretching deformation it remains closed all the time. 

In conclusion, we have presented total energy DFT 
calculations for the stretching deformation of an Al- 
nanowire. The use of a fast local-orbital code has allowed 
us to analyze in detail the nanowire geometrical config- 
urations along the stretching path. At the same time, 
the first-principles tight-binding Hamiltonian obtained at 
each geometry for the electronic system, has allowed us 
to analyze also in a very efficient way, using Keldish- 
Green function methods, the differential conductance of 
the system. Our results reproduce most of the experi- 
mental evidence, and show that the deformation induces 
at particular points sudden rearrangement of atoms that 
manifest itself in jumps of the forces and the conduc- 
tance. We also find that in the final process of deforma- 
tion, the nanowire breaks at a point having a differential 

conductance close to (^TT^j > ana - ^ na ^ the total conduc- 
tance is the result of the contribution of three channels, 
one of them being the predominant one. 
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